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ABSTRACT 

Taking advantage of both the high temporal and spatial resolution of the Atmo- 
spheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), we 
studied a limb coronal shock wave and its associated extreme ultraviolet (EUV) wave 
that occurred on 2010 June 13. Our main findings are (1) the shock wave appeared 
clearly only in the channels centered at 193 A and 211 A as a dome-like enhancement 
propagating ahead of its associated semi-spherical CME bubble; (2) the density com- 
pression of the shock is 1.56 according to radio data and the temperature of the shock 
is around 2.8 MK; (3) the shock wave first appeared at 05:38 UT, 2 minutes after the 
associated flare has started and 1 minute after its associated CME bubble appeared; 
(4) the top of the dome-like shock wave set out from about 1.23 R© and the thickness 
of the shocked layer is ~ 2x10^ km; (5) the speed of the shock wave is consistent with 
a slight decrease from about 600 km s~^ to 550 km s^^; (6) the lateral expansion of 
the shock wave suggests a constant speed around 400 km s~^, which varies at different 
heights and directions. Our findings support the view that the coronal shock wave is 
driven by the CME bubble, and the on-limb EUV wave is consistent with a fast wave 
or at least includes the fast wave component. 

Subject headings: Sun: eruption, Sun: EIT wave, Sun: EUV emission. Sun: CME, Sun: 
shock waves. Sun: radio emission 



Introduction 



In 1997, the EUV Imaging Telesc ope (EIT; iDelaboudiniere et al.l Il995l ) on board the Solar 
and Heliospheric Observatory (SOHO; iDomingo et al.lll995l l observed globally propagating wave- 
like disturbances in the corona, known as EUV waves (also called "EIT" waves), for the first 
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time (jMoses et al.l 119971 : iThompson et alj Il998l ) . Usually, EUV waves are observed as diffuse and 
irregular arcs of increasing coronal emission in the 195 A channel. Compared to Moreton waves, 
E UV waves are a relativ ely frequent phenomenon and their speeds are relatively low. According 
to lBiesecker et al.l (j2002l ). about 7% of the events in t heir large sample displa yed sharp and bright 
wavefronts somewhat reminisc ent of Moreton waves (IThompson et al.l l200d). Such EUV wave s 
are also called "bow waves" by iGopalswamy et al.l (j200Cl ) or "S-waves" bv lBiesecker et al.l ()2002l ). 
In several events t he S-waves coincide s patially with Moreton waves observed at the same time 
( Khan et al.l 120021 : IWarmuth et al.l 12004 ). which would imply that at least S-waves are the long- 
sought coronal counterpart to Moreton waves. It is generally accepted that both Moreton waves 
and S-waves are shock waves. 

However, the physical nature o f the more common diffuse EUV waves is sti l l being debated 



among fast-mode MHD w a ves (e.g.. IWills-Dayev &: ThompsonI Il999l : IWangj l2000l: 



Wu et al 



20011 



Ofman k ThompsonI I2OO2I : ISchmidt k OfmanI l20ld '). slow-mode MHD waves (jWang et al.l 12009) 
and non- wave s related to a current shell or success i ve restructuring of field lines during the c oronal 



mass ejection (jDelannedl2000l : lDelannee et al.ll2008l : IChen et al 



2002 



iEooi; 



Wuelser et al 



Attrill et al. 



200 



3 ). The 



2004:lHoward et al 



observations from the Extreme Ultraviolet Imager (EUVI; 
instr uments on board the twin Solar- Terrestrial Relations Observatory (STER EO, 



20081) spacecraft have advanced our understanding of the diffuse EUV wave (e.g.,lLong et al 



using computational simulation ( Cohen et al.l I2OO' 



2008) 



Kaiser et al 



wave/non-wave hypothesis was fi rst suggested bvlZhukov k Aucherd (120041 ) and has been extended 



I 



Downs et al 



200 



Veronig et al. 


2008 


: ( 


3opalswamv et al. 


2009: 


Attrill et al. 


2009; 


Patsourakos et al 


2OO9I: Ma et al. 


2009; 


Kienreich et al 


2OO9I: 


Zhukov et al. 


2009: 


Veronig et al. 


20ld: Yang k Chen 


20ld). A hybrid 



20 111 ). For detailed discussions 



of the differ e nt models and support in g observations, we ref e r to t he re cent reviews bv IWa.rmuth 
(|200i boinl ): IVrsnak k Clh^ (|2008l ): Iwills-Davev k Attrili (I2OO9I ) and ICallagher k LongI (|20inl )~ 



Due to the scarcity of direct observ ations of shock waves with high spatial resolution, the well- 
know n signatures — type II radio bursts ( Wild k McCreadvlll95C ; I Wild et al.l[l954l :l iNelson k Melrosd 



white light ( Ontiveros k Vourlidasll2009 ) and UV spectra ( Ravmond et 



(Ravmond et a 


.II2OOC: 


Mancuso k Bemporad 


. Klassen et al. 


( 


2000) made a statistical 



2OO9I ) have also been used to study coronal shock waves. iKlassen et al 
study of type II bursts, and found that 90% are associated with an EUV wave. However, the exact 
relationship between the EUV and shock wave is still open because of the low (or even no) spatial 
resolution of the radio data. For the same reason, the origin of the shock wave is also open. Two 
interpretations have been proposed, one suggesting a blast wave ignit ed by the pressure pulse of a 



flare and the other a r guing a piston-driven shock due to a CME (e.g ., iReiner et al.ll200ll : iLin et al 
20061 ; lOh et al.l[2007l : iMagdalenic et ahlboosl ; IVrsnak k Cliveilbood ). 



Recently, IVeronig et al.l (|201Cll ) presented a study of a weak dome-like shock wave observed by 
EUVI which proved the feasibil ity of using EUV observations to study sh ock waves. The Atmo- 
spheric Imaging Assembly (AIA. lTitle k AIA teamll2006l : iLemen et al.ll201ll ) on the Solar Dynamics 
Observatory (SDO) has 7 EUV and 2 UV wavelengths, covering a wide range of temperatures, at 
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high ti me cadence (12s) and resolution (l".4, with 0".6 pixels). Taking advantage of these capa- 
bilities, iLiu et al.l (|20ld ) studied a global coronal EUV wave and found multiple "ripples" for the 
first time; their results support the hybrid EUV wave hypothesis. Direct observations of shock 
waves with high spatial resolution in the low corona would greatly improve our understanding of 
the origin of shock waves and the nature of diffuse EUV waves (hereafter EUV wave) . 



In this paper, we present a detailed analysis of the 2010 June 13 event (cf. iKozarev et al 



including time-dependent ionization and the nature of EUV wave near the solar surface. We 
briefly introduce the data used in this paper in the next section. The observations and results are 
presented in Section 3, and we present our discussion and conclusions in Section 4. 



2. Data and Observation 

On 2010 June 13, a limb solar eruption occurred in the active region NOAA 11079 (S24W91) 
from around 05:30 UT to 05:50 UT. This event involved a filament eruption, an Ml.O flare 
and a CME with a stro ng and short-lived acceleration phase that has been studied in detail by 
Patsourakos et al.l (j2010l ). In particular, a distinct dome-like shock wave associated with this erup- 
tion was observed clearly by the AIA imagers in the coronal channels centered at 193 A and 211 A, 
allowing us to study the morphology and propagation characteristics of the shock wave in detail. 
In addition, a type II burst was also observed by some radio spectrographs, indicating the presence 
of a shock wave. 

To study the shock wave and associated phenomena, we mainly used the following data: 



• The AIA on SDO provides multiple simultaneous high-resolution full-disk images of the corona 
and transition region up t o O.SRf?) above the s olar limb with 1.5 arcsec spatial resolution and 



12 s temporal resolution (jLemen et al.ll201ll ). Seven narrow EUV bandpasses centered on 
specific lines: Fe XVIII (94 A), Fe VIII,XXI (131 A, Fe IX (171 A), Fe XII,XXIV (193 A), 
Fe XIV (211 A), He II (304 A), and Fe XVI (335 A) have been employed. The temperature 
diagn ostics of the EUV emissions cover the range from 6x10^ K to 2x10'' K ( Lemen et al.l 



2011 



Here, we mainly use the channels centered at 171 A, 193 A, 21lA and 335 A (Level 
1.5 images). To show the shock wave more clearly we emplo yed an improved radial filter 
technique first developed by S. Cranmer and A. Engell following iMa et al.l (|2010l ). The radial 
filter technique involves dividing the sun into different concentric "rings" and calculating the 
minimum (Imin) and maximum (Imax) intensity in each radial ring. The scaled intensity at 
each pixel in the ring is then obtained by using the equation lscaied= (I- Imin) /(Imax- Imin). 

The data from EUVI onboard the STEREO Ahead spacecraft (STA). STA is 73.6° ahead 
of AIA on its orbit during the eruption, therefore this event is a limb event for AIA and 
an on-disk event for STA. The two different points of view allow us to construct the 3-D 
structure of the coronal wave and pursue the nature of coronal waves. EUVI observed the 
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chromosphere and low corona in four spectral channels (304 A, 1 71 A, 195 A and 284 A) out 
to 1.7 R0 with a pixel- limited spatial resolution of 1".6 pixel" ^ ( Wuelser et al. 20041 ) . Here 
we only use the 195 A images from EUVI A with a time cadence of 2.5 minutes. 

The data from the solar radio spectrograph of San Vito observatory, which sweeps the fre- 
quency range 25 to 180 MHz every 3 seconds. It monitors solar radio emissions originating 
mainly in the solar corona. It has a low band (25 to 75 MHz) antenna (non-tracking semi- 
bicone) and a high band (75 to 180 MHz) antenna (tracking log-periodic). The radio data 
used in this paper was downloaded from the USAF Radio Solar Telescope Network (RSTN) 
solar radio fixed frequency and spectral data at the NOAA National Geophysical Data Center 



(NGDC) Website ( http : //www . ngdc . noaa . gov ) . 



3. Results 
3.1. Shock Wave in EUV Wavelengths 

3.1.1. Physical Characteristics 

Figure 1 shows the AIA observations (processed using the radial filter technique) of this event 
in the 17lA, 193 A and 21lA channels. The top row images in Figure 1 show the morphology 
of the corona before the shock wave appeared and the bottom row images show the morphology 
during the shock wave propagation. During the eruption, a sharp spherical CME bubble (indicated 
by the black arrows in Figure 1) can been seen in all three channels. For the first time, a spherical 
shock wave (indicated by white arrows in Figure 1) was observed by AIA at EUV wavelengths (see 
the online animations: aia_20100613_171.mpg, aia_20100 613_193.mpg and aia 20100613_211.mpg). 



Note that this event was also which is also studied by iKozarev et al.l (|201ll ). The shock wave 
was ahead of the CME bubble from its first appearance in AIA images at about 1.23 R0 (from 
Sun center) at ~05:38 UT. The thickness of the layer of shocked gas is ~2xl0^ km based on the 
measurement from AIA images along the radial direction (meridional 116 degree) in the 193 A 
channel. 

Interestingly, the shock wave appeared as a distinct bright feature only in the 193 A and 211 
A channels. No clear signature of the shock wave can be identified from the 171 A image in Figure 
1, but the online animation aia_20100613_171.mpg shows that the shock wave appeared as a dark 
feature in the 171 A channel. The 193 A and 211 A channels are dominated by the Fe XII lines 
[log T~6.2] and Fe XIV lines [log T~ 6.3] for AR observations, respecti vely. The 171 A chan nel 



observes the Fe IX line [log T ~ 5.85] for active region (AR) plasma (see lO'Dwyer et al.ll2010l for 
detailed information). We show below that the temperature of the dominant plasma in the shock 
wave is around 2.8 MK. 



The white box in Figure 1 mark the area where we measured the intensity change before and 
during the passing of the shock wave. The intensity tracking in the white box is placed in Figure 2, 
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with the intensities in 211 A and 335 A scaled by multiples 5 and 250, respectively. The plot shows 
that after the shock wave passed, the intensity in the 193 A, 211 A and 335 A channels increased 
rapidly and then went through a rapid decrease followed by a slow one. However, the intensity in 
171 A channel decreased when the shock wave passed and increased slowly until the CME bubble 
arrived. The intensity jumps (maximum intensity divided by the pre-shock intensity) are 1.2, 1.5 
and 1.7 in the 193 A, 211 A and 335 A band, respectively. 



3.1.2. Upward Propagation 

The three dotted lines in the bottom panels in Figure 1 indicate the radial directions at 115°, 
116° and 117 ° clockwise from the solar north pole, respectively, along which we analyzed the 
upward propagation of the shock wave and associated CME. 

Figure 3 contains time-distance slit images showing the shock wave formation and propagation 
in different wavelengths along the three radial directions marked by the dashed black lines in Figure 
1. The images in the left, middle and right columns are taken from the 171 A, 193 A and 211 A 
channels, respectively. For each column, the top, middle and bottom panels are obtained by placing 
the slit along the radial directions at 115°, 116° and 117 °. The dashed curves in each panel of 
Figure 3 are the GOES flux in 1-8 A passband showing the associated flare. A filament associated 
with this eruption rose sharply beginning at 05:32 UT (panel d and g in Figure 3), indicating that 
the filament initiated this solar eruption. 

The CME bubble is most apparent in the 171 A channel, so we use it to study the propagation 
of the CME bubble. In each radial direction, the CME bubble shows similar propagation: going 
through a rapid acceleration phase and then mpging with a nearly constant or slightly decreasing 
speed. The slit image at 117° in the 171 A channel (panel c of Figure 3) is then used to measure 
the distance, speed and acceleration of the CME bubble. The result is given in the top panel of 
Figure 4 with the red, blue and orange lines with plus symbols indicating the distance, third degree 
polynomial fit speed and acceleration of the CME, respectively. The upward speed of the CME 
initially increased from to 500 km s~^ and then decreased to 250 km s~^ before leaving the AIA 
field of view. The corresponding acceleratio n of the CME decreased from around 6 to -4 km s~^. 



(For a detailed analysis of the CME, refer to iPatsourakos et al.ll20ld .) 



The shock wave can be identified in the 193 A and 211 A channels, although its appearance 
differs slightly in each. In the 193 A channel, the shock wave front is relatively sharp and distinctly 
separated from the CME bubble (see panels d-f of Figure 3). However, in the 211 A channel, the 
intensity in the region between the shock wave front and the CME bubble appears evenly increased 
(see panels g-h of Figure 3). As discussed in Section 3.3, this is because it takes a finite time to 
ionize the plasma from Fe XII, which dominates the 193 A channel, to Fe XIV, which dominates 
the 211 A channel. The shock wave first appeared in AIA images at 1.23 R© from the surface of 
the Sun at around 05:38 UT (dotted-dashed line in panel e of Figure 3), two minutes after the 
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CME showed a distinct rise (dotted line) . For measuring convenience, we chose panel e of Figure 

3 to analyze the shock wave propagation along the radial direction. Rough estimates showed that 
the linear fit speed of the shock wave is around 600 km s~^, larger than the CME bubble's linear 
fit speed of 410 km s^^ (not including the rapid rising phase). 

We identified the brightest point of the shock wave in each image and used a semi-automatic 
tracking method to obtain the distance of the shock wave along the radial direction. The result 
is shown in the top panel of Figure 4, with the red, blue and orange lines with diamond symbols 
giving the distance, third degree polynomial fit speed and acceleration of the shock wave front. The 
third polynomial fit speed of the shock wave showed a slight decease from 600 km~^ to 550 km s~^ 
with a monotonically decreasing acceleration less than -1 km s~^. 

The bottom panel of Figure 4 displays GOES soft X-ray flux. The orange curve for 1-8 A, the 
purple curve for 0.5-4 A and the green curve is the time derivative of the 1-8 A flux. The three 
dotted vertical lines in each of panels from left to right in Figure 4 mark the times: 05:36 UT, 
05:37 UT and 05:38 UT, respectively. The Figure shows that the flare began at around 05:35:30 
UT and peaked at 05:39 UT. The distinct disturbance of coronal loops which later became part 
of the CME bubble first appeared at 05:36 UT (left line), while the sharp CME bubble came into 
being at 05:37 UT (middle line). The shock wave appeared at 05:38 UT (right), which is later than 
the first appearance of both flare and CME. 



3.1.3. Lateral Propagation and EUV Wave 

Figure 5 gives base difference images showing the evolution of the shock wave and the associated 
EUV wave. The EUV wave here refers to the disturbance propagating laterally along the solar 
surface which normally appears as a projected bright circular feature propagating on the solar disk 
in EUV difference images. Panels a-h of Figure 5 are AIA base difference images in 193 A channel 
obtained by subtracting the 193 A image at 05:32:08 UT from the present images. Panels i-1 of 
Figure 5 are EUVI A images in the 195 A channel obtained by subtracting the 195 A image at 
05:33:00 UT from the current images. The AIA data (panels a-h, see also the online animation 
basedif_aia_20100613.mpg) show a close relationship between the shock wave and the EUV wave: 1) 
they appeared at the same time ~ 05:38:08 UT; 2) the front of the laterally expanding shock wave 
is tightly connected to the position of the EUV wave front. In EUVI observations, we could not 
identify the coronal shock wave, probably because the shock wave is weak and the column depth 
is small. The EUV wave can be identified in EUVI as a circular disturbance propagating from AR 
11079 (see also the online animation EUVI_A_20100613.mpg). 

The top panel of Figure 6 is a projected base difference image (see the online animation 
project_diff_20100613_193.mpg for more information) where the horizontal coordinate is polar angle 
(clockwise from north pole of the Sun) and the vertical coordinate is radial distance in solar radii. 
The five horizontal lines (which are parallel to the solar surface) from bottom to top indicate the 
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layers located at 0.98, 1.01, 1.04, 1.07 and 1.11 R©, respectively. By placing the slit along these 
layers, we obtained five different slit images from the base difference images in 193 A. These slit 
images are placed in panels b-f of Figure 6. The horizontal coordinate of the slit images is the 
meridional angle in degrees while vertical coordinate is time in minutes. The black area indicates 
the dimming area changing with time and the white features indicate the EUV brightening. Taking 
the propagating white feature as the coronal wave front, we estimate the EUV wave speed. The 
dotted lines and the numbers around them in Figure 6 indicate the linear fit distances of the EUV 
wave fronts and the estimated speeds. The propagation of the EUV wave towards the north (left) 
and south (right) is not symmetric; the speed in the southward direction appears slower than the 
northward direction. 

A bright feature was also observed propagating in the reverse direction, implying the reflection 
of the EUV wave (indicated by the "RF" in panel d and e of Figure 6). However, a detailed analysis 
of this observation is left for future work. 



3.2. Shock Wave in Radio Observation 

According to the National Geophysical Data Centers event listing, a type II burst associated 
with this event was observed by several radio spectrographs. In Table 1 we list the extracted dates, 
observatory stations, event start time, event end time, spectral class (SC), lower frequency (LF), 
upper frequency (UF) and estimated shock speeds in km s~^. The start time of the type II fits well 
with the AIA observation of the shock wave and there is no other solar eruption at the same time. 
Therefore, the shock wave observed at radio and EUV wavelengths should be the same one. 



Table 1: 


Record of the type II occurred 


on June 13, 


2010. 




Station 


Start time 


End time 


SC 


LF(MHz) 


UF(MHz) 


Shock Speed (km/s) 


BLEN 


0537.0 


0541.1 


II 


175X 


332 




CULG 


0537.0 


0543.0 


II 


57X 


150 




HIRA 


0537.0 


0550.0 


II 


50 


310 




CULG 


0538.0 


0551.0 


II 


57X 


200 


700 


SVTO 


0538.0 


0552.0 


II 


35 


180 


621 


LEAR 


0539.0 


0548.0 


II 


55 


180 


665 



Figure 7 shows the metric type II radio emissions from San Vito radio (SVTO) spectrograph. 
The top panel is the dynamic spectrogram. The type II burst occurring from 05:38 UT to 05:53 UT 
can be easy identified. The black curve indicates the harmonic frequency (2/p) in the type II burst 
along which we measured the local plasma frequency (fp) drift. The result showed that 2fp drifts 



4tp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_RADIO/SPECTRAL/2010/SPEC_NEW.10 
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from 165 MHz to 40 MHz. Considering the relationship between plasma frequency and density 

fp = 8.98 X 10^^/^, (1) 

we determine the density of the shock wave front (middle panel of Figure 7). The density of the 
shock wave dropped from 8.6 x 10'' cm~^ at around 05:38 UT to 4.5 x lO^cm"^ at about 05:52 UT. 



The i nset in the middle panel is the coronal plasma density model from ISittler &: Guhathakurta 
(|l999l ) as used bv lLin et al.l tood ). 



n{z) 



no ai z'^ ^ [1 + 03 z + 04 + as z^] 



(2) 



z = l/{l + y),ai = 0.001272, 02 = 4.8039, 

as = 0.29696,04 = -7.1743, ag = 12.321 

where y is the height above the solar surface in solar radii and no is the electron number density at 
the solar surface, which is chosen as 10^ cm~^ here. The density model was used to estimate the 
position of shock wave front (the black line in the bottom of Figure 7), showing that the shock wave 
formed at around 1.25 Rq. The derived speeds are displayed by the dotted and dashed lines in the 
bottom panel of Figure 7. The dotted line shows the linear fit speed (527 km s~^) and the dashed 
line shows the second degree polynomial fit speed which decreased from around 600 km s~^ to 400 



km s~^. The speed of the shock wave estimated using the model of Sittler Sz Guhathakurta (jl999l ) 
is smaller than the speeds listed in Table 1. As the speed derived from radio emissions is strongly 
dependent on the coronal density model used, the error of the speed is quite large. Considering 
this situation, the speed of the shock wave measured from radio emission is consistent with that 
obtained from AIA observation. 



The two black stars mark the frequencies used to measure the density jump at the shock using 
the equation, 



(3) 



As described in lVrsnak et al.l (j2002l ). the band-split frequencies in type II emission map the electron 
densities behind and ahead of the shock front. In front of the shock the plasma is characterized by 
the electron density ni and emits radio waves at the frequency fi (lower frequency branch) while 
the plasma behind the shock is characterized by the electron density n2 and emits radio waves at 
the frequency fjj. At 05:40 UT, the lower and upper frequencies of the harmonic bands in type 
II burst are around 132ib5 MHz and 165ibl5 MHz (marked by the black stars in the top panel 
of Figure 7), respectively. U s ing eq uation (3), we obtained the density jump X = 1.56 ± 0.1 at 
05:40 UT. iGopalswamy et al.l (120111 ) obtained a similar result using radio data from HiRAS. For an 
oblique shock (taking the adiabatic index 7 = 5/3) the Alfven Mach number Ma and the density 
jump X are related by: 



(M| - Xf[5f3X + 2M\ cos^ e{X - 4)] 
+M\ X sin^ 0[(5 + X)M\ + 2X{X - 4)] = 



(4) 
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(jVrsnak et al.lbood 'l. For /? ^ 0, in the case of the perpendicular shock (6* = 90°), MA=y/X{X + 5 + 5/3)/2(4 
and in the case of longitudinal shock (9=0°), M^=\/X=1.25. Considering Ma=v/va and Ma= 



1.35±0.1 (the actual error maybe larger) this produces an estimated Alfven speed 450ib30 km s~^. 
The estimated sound speed Cg = {jKT)/m is 156ib30 km s~^ for a pre -shock temperature o f 
(1.8±0.4) X 10^ K. According to the DEM solutions for regions 2 and 3 in iKozarev et al.l (|201ll ). 
the pre-shock temperature (peak value) is around logT = 6.25 it 0.1 (1.8ifc 0.4 MK). The a rea fo - 



X)- 



cused on in this paper differs slightly from the region 2 and 3 studied by iKozarev et al.l ()201. 



so there may be a small discrepancy in the temperature, although it should not be too different. 
In the following context, we using 1.8ib0.4 MK as the pre-shock temperature T^. The resulting 

-1 



fast magnetosonic speed Cf 



+ c2 is 476±38 km s" 



In the solar atmospher e with /2=0.6 



and 7=5/3, the magnetic field strength can be derived (from the equation 2.48b of |Priestlll982l ) 
n^J"^ Va- For a pre-shock density of 6x1 0^ cm~^ from the type II frequency 



as B = 3.57 x 10" 



at around 05:40 UT, the value of va implies B ~ 1.3 Gauss. iGopalswamv et al.l (|201lh obtained 
an estimate for B that is consistent with this using a completely different technique based on the 
standoff distance between the flux rope and the shock. 



3.3. Verifying the Shock Wave Interpretation 



Coronal shocks have been identified from white light coronagraph imag es by morphology and 



density enhancement ( e.g. 



mavmond et al 



200C 



Vourlidas et al 



2OO3I: lOntiveros k. Vourlidasll2009l ) and from UV spectra 



Mancuso &: BemporadI |2009| ) by measuring oxygen kinetic temperatures. 



AIA shows brightening ahead of the bright CME loops, but that could in principle be part of 
ejected streamer rather than a shock. The typical attributes of a shock are density compression, 
gas heating and subsequent gas ionization. The intensities obtained from AIA imagers combine 
the electron density and ionization state of the plasma. To test the relation between compression 
and heating by an MHD shock, we employ Vs=600 km s~^ from observed speed of shocked front 
and /02//9i=1.56 from the band splitting in the radio. If the shock was purely gas dynamic, with 
a sound speed Cs=156 km s~^, and a Mach number M=600/156=3.8, the compression would be 
P2//Oi=3.3, instead of the observed 1.56. Therefore the magnetic fiel d B is important. Acco rding 
to the jump condition for a perpendicular shock (Equation 2.19 of iDraine &: McKed Il993l ). the 
compression ratio is 

P2^ 2(7 + 1) 

Pi D + [Z)2 + 4(7-M)(2-7)M/]V2' 
where L> = (7 - 1) + (2M~2 ^ ^MJ^). Assuming M=3.8 and /)2/pi=1.56, we find Ma=1.55. 



(5) 



For a perpendicular shock, the shocked plasma parameters (^2, P2, P2, B2) are related to those 
of the unshocked plasma (fi, pi, pi, B\) by the equations for conservation of momentum and 
energy, 

P2 + Bl/{2ii) + p2vl =pi + Bl/{2p) + pivl (6) 

,2 



(|Priesdll982l 1. Considering ps/pi = X, B2/BI = X, v^jvx = l/X, p 



nkT, 



7pl//jl, 
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v\ = M = vl/cs and Ma = v1/va, the relationship between the pre-shock temperature 

and post-shock temperature can be obtained as 

^2 = l[l + (1 _ 1 _ ^!^)^m2]. (7) 

Taking X=1.56, M=3.8, Ma=1.55 and 7 = 5/3, T2/T1 is around 1.57 according equation (7). 
Assuming the pre-shock temperature of Ti=(1.8±0.4) x 10^ K as we explained in Section 3.2, the 
post-shock temperature T2 is around (2.8=b0.6)x 10^ K. 

Using the ionization rates qi at T=2.8xl0*' K, we can estimate the ionization time scales and 
then compare them with observed time scales of the band ratios. The ionization time scale to reach 
Fe XIV (211 A) and Fe XVI (335 A) can be obtained from the following equations: 

t2ii = + ) (8) 

QFeXII QFeXIII 

1.1 1 lis 
i335 = — ( + + + )• (9) 

ne qPeXII qPeXIII qPeXIV QFeV 

Where qpeXii, QFeXiii, QFeXiv and qpexv are ionization rate coefficients in units of cm^ s~^. Given 
the temperature of plasma the ionization rate coefficients can be easily obtained from CHIANTI by 
using the code ioniz_rate.pro. In Table 2 we list the ionization rate coefficients at 2.4 MK, 2.6 MK, 
2.8 MK, 3.0 MK, 3.2 MK, respectively, in unit of cm''' s'^. Taking ne=6x 10^x1.56=9.4x10'^ cm"^, 
the ionization time scales can be obtained by using the equations (8) and (9). The corresponding 
results are placed in the right two columns in Table 2. 

Table 2: The ionization rate coefficients (in unit of cm'"^ s^"*^) at different given temperatures and 
the corresponding estimated ionization time scales. 



T (MK) 


ClFeXII 


QFeXIII 


QFeXIV 


QFeXV 


t21l(s) 


t335(s) 


2.4 


4.0982035e-10 


2.4158748e-10 


1.4844586e-10 


1.0106319e-10 


110 


386 


2.6 


4.6829123e-10 


2.8404318e-10 


1.7771562e-10 


1.2294273C-10 


94 


323 


2.8 


5.2553271e-10 


3.2686128e-10 


2.0770653e-10 


1.4567346e-10 


83 


277 


3.0 


5.8123741e-10 


3.6961965e-10 


2.3806008e-10 


1.6894433e-10 


74 


242 


3.2 


6.3520065e-10 


4.1199241e-10 


2.6848624e-10 


1.9249588e-10 


66 


215 



Figure 8 displays time intensity ratio profiles, which are derived from the intensity changes 
in different wavelengths (Figure 2) for the region marked by the white box in Figure 1. The rise 
phase of the intensity ratio profile is indicating the ionization time. The ionization time scales 
of 211 A and 335 A obtained from observation are around 100ibl2 seconds (marked by the top 
two vertical lines in Figure 8) and 275±12 seconds (marked by the bottom two vertical lines in 
Figure 8), respectively. In general, the ionization times obtained from observation are consistent 
with that derived from the theory of shock wave, 85ibl3 seconds for t2ii and 288ib53 seconds for 
t335. In other words, the brightening ahead of the CME bubble is consistent with the shock jump 
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conditions, while it would be purely coincidental if the structure were not a shock. The theoretical 

estimates of the ionization times arc limited by the uncertainty in the pre-shock temperature, and 
this might be improved by further analysis of the AIA data. We also note that the intensity ratios 
of the raw band intensities can be misleading, because they convolve the ionization times with 
the time scale for the increase of emission from the shocked gas compared to the intensity from 
the foreground and background plasma outside the shock. We have estimated the foreground and 
background emission based on the variation of the 171 A band, and in this case the time scales 
for the intensity ratios are similar to those obtained from the raw band intensities. Moreover, a 
more detailed analysis would include the pre-CME density profile and the adiabatic cooling as the 
shocked gas expands. It should be possible to obtain accurate values of the electron temperature 
at different positions around the shock to study electron-ion temperature equilibration and energy 
losses to particle acceleration. 

4. Discussion and Conclusions 

On June 13 2010, a CME-associated shock wave was directly observed by AIA. Combining 
these observations, EUV data from STEREO A and radio data from SRS, we present a detailed 
study of the shock wave and its low coronal reactions. 

4.1. Main Findings 

Our main findings are as follows: 

• The shock wave front first appeared at 05:38 UT as a dome-like bright feature, which is 
clearly distinguishable from its associated CME bubble. The top of the dome-like shock wave 
originated from about 1.23 Kq with the thickness of the shocked layer around 2x10^ km. A 
type II radio burst occurred at the same time, further confirming the presence of the shock 
wave. The electron density jump X (the ratio of the plasma density behind the shock and 
that in front of the shock) is around 1.56 indicating that the shock wave is a weak shock 
wave. The pre-shock density and the density jump from the radio data, along with a shock 
speed of 600 km s~^, imply a post-shock temperature of 2.8 MK. The ionization time scales 
of 211 A and 335 A obtained from AIA observation are consistent with the ionization time 
scale derived from shock wave jump relation in general, which proves the dome-like bright 
feature is a shock wave. 

• The shock wave appeared distinctly only in the channels centered at 193 A and 211 A with 

the estimated relative intensity of the shock wave front increasing by 1.2, 1.5 and 1.7 in 193 
A, 211 A and 335 A, respectively. All of the EUVI channels showed an increasing intensity 
with the exception of the 171 A channel which showed a clear intensity decrease. 
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• The upward speed of the shock wave shows a shght decrease from about 600 km s~ to 550 
km with a deceleration less than 1 km s~^. The lateral speed of the shock wave shows no 
acceleration, but varies according to both height and direction. The linear fit lateral speeds 
towards the north range from 246 km s~^ to 397 km s~^ while the speeds towards the south 
range from 342 km to 486 km s^^. The lateral speed of the shock wave seems to increase 
with the height below 1.1 Rq. In general, the lateral speed of the shock wave is less than its 
upward speed. 



4.2. Origin of the Shock Wave 

The high time cadence of AIA observations makes the pursuit of the origin of the low coronal 
shock waves possible. AIA data shows that a filament rising was observed at 05:32 UT, at the very 
beginning of this solar eruption. About four minutes after the filament rising, a flare and CME 
bubble appeared nearly simultaneously. At 05:38 UT, two minutes later after the flare and CME 
bubble appeared, the shocked front came into being. Considering the time sequence, it is likely 
that the filament is a trigger for the whole eruption. Because the fiare and CME bubble appeared 
simultaneously and earlier than the shock wave, neither of them can be eliminated as the driver of 
the shock wave according to the time sequence alone. However, considering the strong similarity of 
the shock wave and the CME in both morphology (dome- like) and kinematics (slightly decreasing), 
we would suggest that the CME bubble played the role of a piston in driving the shock wave. 



4.3. Nature of the EUV Wave 



The direct observation of the low coronal shock wave by AIA strongly supports the existence 
of a dome-like fast mode wave observed in EUV wavelengths. The fastest front of the EUV wave 
coincides with lateral expansion of the shock wave. Therefore our result supports the hypothesis 
that the EUV wave is a fast wave or at least includes a f ast wave c o mpon ent. Considering the pro- 
jection effect of the on-disk observation as discussed in iMa et al.l (|2009l ). the on-disk observation 
may include non-wave component. However due to the weakness of on-disk wave front obse r ved b y 
EUVI, we did not get strong evid ence for non-wave component as proposed by iLiu et al.l (j20ld ). 
More recently, IChen Wul (j201ll ) showed an on-disk event and found both fast and slow propa- 
gation fronts. However, there are still inconsistencies between the observed phenomena and the 
different theories proposed to explain the EUV wave. The high cadence observations from SDO 
should allow these issues to be resolved. 
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Fig. 1. — AIA images showing the morphology of the shock wave (reversed color table). The AIA 
images were contrast enhanced using a radial filter technique described in Section 2. (Animations 
aia_20100613_171.mpg, aia_20100613_193.mpg and aia_20100613_211.mpg are available in the online 
journal.) 
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Fig. 2. — Intensity flux tracking in the white box in Figure 1. The green, pink, blue and ohve 
curves refer to the intensity flux in 171 A, 193 A, 211 A (x5) and 335 A (x250). The dashed black 
line and dash-dot line refer to the time when the shock and CME bubble arrived in the white box. 
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Fig. 3. — AIA slit images (reversed color table) in 171 A (left column), 193 A (middle column) and 
211 A (right column) showing the time-distances of the CME bubble and shock wave along 115° 
(top row), 116° (middle row) and 117° (bottom row), respectively. The contoured dash curves are 
the GOES 1-8 A flux. 
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Fig. 4. — CME dynamics and GOES flux. The red, blue and orange curves with plus (diamond) 
symbols indicate the time-distance, time-speed, and time-acceleration profiles of CME bubble 
(shock). 
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Fig. 5.— AIA 193 A (top two rows) and EUVI A 195 A (bottom row) observations show 
the propagation of shock wave and EUV wave. Animations (basedif_aia_20100613_193.mpg and 
EUVI_A_20100613.mpg) are available on line. 
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Fig. 6.— AIA 193 A slit images showing the coronal wave propagating at different heights 
parallel to the solar surface. The base time of the event is 05:32:08 UT. Animation 
(project_diff_20100613_193.mpg) is available on line. 
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Fig. 7. — The dynamic spectrogram from San Vito observatory in the range 25-180 MHz and the 
derived density and speed of the shock wave. 



-24- 




Fig. 8. — Time intensity ratio profiles from the white box in Figure 1. The thick black curve for 
the intensity ratio of 211/193 and the thin black curve for the intensity ratio of 335/193 (x30). 



